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Abstract 

The Jeans scale is estimated during the coexistence epoch of quark-gluon and 
hadron phases in the first-order QCD phase transition. It is shown that, contrary 
to previous claims, reduction of the sound speed is so little that the phase transition 
does not affect evolution of cosmological density fluctuations appreciably. 
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As is well known, the Jeans length scale, above which density fluctuations can grow, 
is given by the sound speed Cg multiplied by the dynamical time scale. In the radiation- 
dominant stage of the early universe, Cg is equal to 3~^/^ times the light speed, which we 
take to be unity here, so density fluctuations can grow only on super-horizon scales. 

Recently, however, Jedamzik [Q] and, independently, Schmid, Schwarz and Widerin 
[0] proposed that if a cosmological quark-hadron phase transition is of flrst-order, which 
is supported by some of the recent lattice simulations 0], the sound speed effectively 
vanishes during the coexistence period of the quark-gluon and the hadron phases. If this 
were true, the Jeans scale would become vanishingly small, and even sub-horizon-scale 
fluctuations should grow during this epoch. This temporal decline of the Jeans length 
would imprint a nontrivial feature on the processed spectrum of fluctuations around the 
solar-mass scale. @] Furthermore, it has been suggested that such anomalous growth of 
perturbations could induce efficient formation of primordial black holes (PBH's)[0]. The 
resultant mass function of the PBH's would have a peak on the solar-mass scale which is 
about the scale of MACHO's. Thus the idea raised in Refs. 1) and 2) is very interesting 
and may be astrophysically important. Unfortunately, however, the actual evolution of 
the sound speed, as well as that of the Jeans length, has not been given in either paper, 
and their claims are limited to a qualitative level. 

In the present paper we present the result of an explicit calculation of the sound speed 
using the bag model We first review how the quark-hadron phase transition proceeds 
in this model and then calculate the sound speed under two distinct conditions. 

In the bag model |] the energy density of the quark-gluon phase, pg, is the sum 
of that of the relativistic particles and the QCD vacuum energy density, namely, the bag 
constant, B. 

p,(r) = ^^,r^ + i?, (1) 

where Qq is the effective degree of relativistic freedom. Qq = 51.25 when the number of 
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species of relativistic quarks is two, and Qq = 61.75 when it is three. The pressure, Pg, 
and the entropy density, Sg, are given by 

P,{T) = ^9,T'-B, (2) 
s,iT) = ^,,T^ (3) 

respectively. On the other hand, the energy density, p^, the pressure, ph, and the entropy 
density, Sh, of the hadron phase are expressed as 

Ph{T) = "^g.T', (4) 
Ph{T) = ^9hT\ (5) 
sdT) = ^9hT\ (6) 

where Qh — 17.25 is the degree of relativistic freedom in this phase. At the coexistence 
temperature, Tc, we have Pq{Tc) = ph(Tc), and hence it is related to the bag constant by 

B = ^{9,-9h)T^. (7) 

In this model, the phase transition from the high-density quark-gluon phase to the 
low-density hadron phase starts with bubble nucleation, whose probability per unit time 
per unit volume, P(T), is given by 

where L is the latent heat per unit volume of the phase transition, a is the free energy 
per unit surface area of the bubble, C is a constant, and r], which represents the degree 
of super cooling, is defined by 

T — T 

The exponent of P(T) is singular at T^, and P(T) vanishes then, but it increases drastically 
when T drops slightly from Tc. 



Hence the phase transition is expected to proceed as follows. In the course of the 
cosmic evolution the cosmic temperature decreases, and until T = the entire universe 
is occupied by the quark-gluon phase. At this temperature, the hadron phase does not 
appear yet, and the super cooling (T < T^) occurs. But only a small degree of super 
cooling, 1] ~ 10^'^, actually takes place since P{T) soon becomes large^]. Then the 
nucleation of the hadron phase starts, and the universe is reheated by the latent heat 
up to Tc so that the pressure equilibrium between two phases is realized and subsequent 
nucleation is suppressed. Afterwards, the phase transition proceeds through expansion of 
nucleated bubbles whose speed is small and estimated to be around 1O~^0, keeping the 
balance between heating due to continuous liberation of latent heat and cooling caused 
by cosmic expansion. The temperature and the pressure remain constant in this stage, 
but the mean energy density. 



decreases gradually from PqiTc) to Ph(Tc) as the volume fraction of the quark-gluon phase, 
fg, drops. During this epoch the scale factor expands by a factor of about {gq/ghY^^ — 1-4. 

It is in this coexistence regime that the authors of Refs. 1) and 2) suggested the sound 
speed vanishes. According to Ref. 1), an adiabatic compression of the system in the 
equilibrium two-phase mixture induces conversion of the low-energy hadron phase into 
high-energy quark phase, and then the mean energy density grows due to the change of 
fq, with the pressure remaining constant. This implies that pressure response to a change 
in energy density is anomalously small, and hence, by their argument, the sound speed 
should effectively vanish. 

In order to clarify whether such an intuitive discussion is correct, we must first find a 
formula for the sound speed in an inhomogeneous medium consisting of the mixture of two 
different phases. While such a problem has not been fully discussed in any cosmological 
context, it is a rather familiar issue in the field of engineering, e.g., propagation of sound 




(10) 
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waves in a water- vapor system in a boiler etc[]^. 

In the rest frame of the sound wave front]], suppose that the medium ahead of the 
wave front is characterized by a pressure p, the temperature T, energy densities pq and 
Ph, volume fraction of the quark-gluon phase fq, and flow velocity Cg, and that the medium 
behind it is characterized hj p + dp, T + dT, pq + dpq, pu + dph, fq + dfq, and Cg + du, 
respectively. Then the continuity equation reads 

fq (P + Pq) l^Cs + {1- fq){p + Ph) 7^C, 
= {fq + dfq) {P + Pq + dp + dpq) 7^ (Cg - du) 

+ (1 - /q - dfq) {p + ph + dp + dph) 7^ {cs - du) , (11) 

per unit area, where 7 = (1 — c^)~^/^ and 7 = {1 — (c^ — du)'^} are the Lorentz factor 
0. The momentum equation is given by 

P+fq{P + pq) I^CI + (1 - /,) {p + Ph) I'cl 
= p+dp+{fq + dfq) {p + Pq+dp + dpq) 7^ (c^ - du)"^ 

+ (1 - /, - dfq) {p + ph + dp + dph) f (c, - duf . (12) 

From (0) we find 

^ 1 - eg /g {dp + Pq) + (1 - fq) (dp + dph) + (pq - ph) dfq ^ 1 - 0^ dp + dp 

""""^I + Cg fqiP+Pq) + il-fq)iP + Pk) " 1 + P + p ' 

(13) 

Inserting this into (0), we obtain 




t Since the time scale of QCD is much shorter than the cosmic expansion scale, the following discussion 
without the effect of cosmic expansion applies to our problem unchanged. 

■''We thank Dr. Schwarz for pointing out that the nonrelativistic equations were employed in the 
original version of our paper. 
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Thus the sound speed in a mixture is obtained by simply replacing the energy density p 
with the mean value p in the usual formula for a single-component system, 



dp, 



(15) 



The pressure derivative of p is written by 



dp dpn dph , . dfq 

d^ = ^^^ + ^^-^^^^ + ^P^-P'^d^- 

When the sound speed in the hadron phase and the quark-gluon phase equals to that in 

the pure radiation, the sum of the first and second terms of the equation (ITB]) provides 



the value of 3^^/^. Hence the difference from the single- component system is given by the 
third term and since {pq — ph) is positive, if the volume ratio of the quark phase increase 
as the pressure becomes higher, the sound speed in the complex phase should be reduced. 
In order to evaluate the degree of reduction, the estimation of dfq/dp is necessary. 

The next question is under what conditions we should calculate the derivative (|T^). 
Usually one evaluates the sound speed under adiabatic conditions, because a change of 
the entropy associated with sound-wave propagation is small in most cases of interest. 
Since this is also the condition Jedamzik proposed to adopt [|l|, let us first calculate the 
sound speed under this condition. 

The isentropic condition reads 

fqSg-fCs + (1 - fq) Sh^Cs = {fq + dfg) {Sg + dSq) 7 (c, - du) 

+ fg- dfg) {sh + dsh) 7 {cs - du) . (17) 
Then using the average entropy density, s, defined by 

S = fqSq + (1 - /,) Sh , (18) 

and combining this with the continuity equation (|l3)j we find 

dp ds 



p + p s 
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(19) 



Using the relation p = sT — p, we can see this is nothing but the second law of thermo- 
dynamics which gives us no information about dfq/ dp. The adiabatic condition is useless 
for the purpose of calculating the sound speed qualitatively. Hence we must find a more 
appropriate one, understanding the dynamics of the phase transition better. 

What is often adopted in the literature apart from the isentropic condition is the 
conservation of the quality parameter, x, which is the energy fraction of the gas component. 
The theoretical sound speed in a water- vapor system evaluated under this condition agrees 
with experimental values qualitatively!]^. In the present case the corresponding quantity 
may be defined by the energy fraction of the high-energy quark-gluon phase, Xq, as 

= , Jn\, ' (20) 

fqPq + [^- fq)ph 

with which dfq/ dp is expressed as 

n._ f (^ - f \ rln f (^ - f \ rlT 

(21) 



dfg fq (1 - fq) dph fq (1 - fq) dpq fq (1 - fq) dXq 



dp Ph dp Pq dp Xq{l—Xq) dp 

Since we cannot measure the sound speed experimentally in our case unlike in a water- 
vapor system, we must work out the appropriate value of dxq/dp from a theoretical view 
point alone. Here we propose to calculate the sound speed under the condition dxq/dp = 0. 
This is appropriate because the transition of the phases through bubble nucleation is 
totally suppressed at the coexistence temperature, as seen in (H), and the expansion 
speed of bubbles is so small that energy transfer through bubble expansion or contraction 
is also expected to be negligible during sound-wave propagation. 



Therefore we substitute this condition into ([Tq) and find 



dp 
dp 



3p(k + LJl] . (22) 

\Pq ph J 



We thus obtain the sound speed under the constant-quality condition as 



71 r?-^ 



w 1 



fq--.] +4 



1 . (23) 



Here y is defined by 

..^i^^tol. (24) 

Ph[ic) igh 

where we have used (^. Since y + 1/y > 2, we find from (^31) that the sound speed 
in the mixed state is indeed smaller than that in the case of pure radiation and it can 
be arbitrarily small if y is much larger or much smaller than unity. In the present case, 
however, the actual value of y is given hj y = 3.63 for = 51.25 and y = 4.44 for 
Qq = 61.75. Hence reduction of the sound speed is mild, and even the minimum sound 
speed, which is realized at fq = 0.5, is as large as 

cf " = (0.77 ~ 0.82) X ^ (25) 

V3 

for Qq = 61.75 and 51.25, respectively. 

Thus the intuitive claim that ~ is proved to be inappropriate as long as we consider 
the sound speed at the coexistence temperature. Although it is true that the coexistence of 
different phases at the first order phase transition reduces the velocity of the sound wave, 
its degree in the case of the quark-gluon to hadron QCD phase transition is insignificant 
and has no drastic effect on the development of cosmological density perturbations. 

Turing back to the original expression ([l^), we see that the only way to realize a 
vanishingly small sound speed is to have an extremely large value of dfq/dp. But this is 
not possible at the coexistence temperature since the nucleation rate (^ vanishes then. 
Nevertheless it may be possible to realize a large value of dfq/dp at a smaller temperature 
when the nucleation rate rises significantly, and the sound speed may become vanishingly 
small then. Unfortunately, however, in the actual evolution of the universe, in the bag 
model the degree of super cooling is expected to be very small, as mentioned above, so 
even if the sound speed vanished in this era, the duration of this era is too short to leave 
any trace in the spectrum of density fiuctuations. 

Until now we have assumed that the QCD phase transition is of first order and subcrit- 
ical fiuctuations are negligible. If, on the contrary, it were of second order, and subcritical 
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fluctuations were effective, the phase transition would proceed while maintaining an equi- 
librium population between quark-gluon and hadron phases. Under such a circumstance, 
fq would change only mildly against perturbation by a sound wave. Thus dfq/dp cannot 
take a large enough value to suppress the sound speed drastically in these cases. 

We therefore conclude that the quark-hadron phase transition has no appreciable effect 
on the final spectrum of density fluctuations. 
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